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ABSTRACT 

We study the kinematics of GALEX-selected Ha knots in the outer disk (beyond R25) of NGC 
628 (M74), a galaxy representative of large, undisturbed, extended UV (Type 1 XUV) disks. Our 
spectroscopic target sample of 235 of the bluest UV knots surrounding NGC 628 yielded 15 Ha 
detections (6%), roughly the number expected given the different mean ages of the two populations. 
The measured vertical velocity dispersion of the Ha knots between 1 — 1.8i?25 (13.5 — 23.2 kpc) is 
< 11 km s~^. We assume that the Ha knots trace an 'intermediate' vertical mass density distribution 
(between the isothermal sech^(2) and exponential distributions) with a constant scaleheight across the 
outer disk {hz = 700 pc) and estimate a total surface mass density of 7.5 Mq pc~^. This surface mass 
density can be accounted for by the observed gas and stars in the outer disk (little or no dark matter 
in the disk is required). The vertical velocity dispersion of the outer disk Ha knots nearly matches 
that measured from older planetary nebulae near the outskirts of the optical disk by Herrmann et al., 
suggesting a low level of scattering in the outer disk. A dynamically cold stellar component extending 
nearly twice as far as the traditional optical disk poses interesting constraints on the accretion history 
of the galaxy. 

Subject headings: galaxies: individual (NGC 628) - galaxies: star clusters - galaxies: structure 



1. INTRODUCTION 

The outskirts of a galaxy are expected to host sig- 
natures of disk formation and hierarchical accretion 
because of the l ong dynamical times at these radii. 
iToth fc Ostrikerl (fl992;l argued that the thinness and 
coldness of inner galactic disks (inside the optical radius, 
^25) posed significant problems for hierarchical models of 
galaxy formation. Qualitatively, these concerns become 
more pronounced if one can establish that cold galactic 
disks extend to even larger radii. 

Various studies have revisited the cold disk prob- 
lem, generally finding that accretion events are less 
destructive than o r iginall y envisioned. Most recently, 
iKazantzidis et ahl ()2009 l) studied the dynamical re- 
sponse of thin galactic disks to bombardment by cold 
dark matter substructure out to large radii in fully self- 
consistent, dissipationless iV-body simulations. They 
found that disks survive these bombardments, but that 
they produce considerable thickening and heating at 
all radii, substantial flaring, and an increase in the 
stellar surface density in the disk outskirts (the lat- 
ter due to outward radial migration of old stars dur- 
ing the growth and redistribution of disk angular mo- 
mentum). O bservations of the dynamical state of outer 
stellar disks (jChristlein fc Zaritskv I l2008f ) demonstrate 
that outer disks generally continue to obey the flat ro- 
tation curves of inner disks, with no increase in the in- 
plane velocity dispersion. Here we present the face-on 
kinematics of one nearby galaxy, NGC 628 (M74), and 

* This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 



compare to both neutral hydrogen and existing stellar 
measurements at smaller radii. NGC 628 is a prototypi- 
cal Grand Design spiral galaxy (type SA(s)c, dynamical 
mass within the studied region ^ 3. 3 x 10^^ Mp), assum- 
ing vrot = 200 km s~_[ out to 2.3i?9.R: lThilker et al. \\200% 
iKamphuis fc Briggs II199 2V and is by far the largest and 
most massive member of its small group (the brightest 
member after NGC 628, UGC1176, is - 4.5 mag fainter 
and over 125 kpc away). NGC 628 shows a standard ex- 
ponential optical light profile to -R25, with only a very 
slight possible downbending in the profile to ^ 1.3i?25 
the e xtent of the deep optical observations (iNatali et al 



Using ultraviolet (UV) imaging, iThilker et al 



classified NGC 628 as a Type 1 extended UV 
(XUV) disk, due to the structured, UV-bright emission 
complexes seen in the outer disk. Because the stellar 
disk appears largely undisturbed (both in optical and 
UV imaging, and, as we will show, from the kinematics) 
and because NGC 628 dominates the dynamics of its lo- 
cal environment, we treat it as representative of large, 
isolated spiral galaxies in the nearby universe. 

The outer disks of galaxies have r eceived much 
recen t attention, both obser v ational (' Thilker et al. 
20071 : iGil de Paz et al. I |200 ,7l: iZ aritskv fc ChristleiiT 
2007; Christlein fc Zaritskv 2008; Herbert-Fort et al 



2009; .Truiillo et al. , .2009; .Herrmann et al. I I2OO90 an d 



theoretical (iBush et al.j 120081; iR oskar et al. I |2008a_, 
IKazantzidis et al. Il2009[ ). The surge in interest in outer 
disks has been fueled by recent ultraviolet (UV) ob- 
servations of nearby disks with the GALEX satellite 



servations oi nearoy 
(IMartin et al. 1 120051; 



iFerguson et al 



Thilker et a"ni2007D . Previously, 



( 1998f ) had used deep Ha imaging to dis- 
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cover star formation in an extended component around 
three nearby galaxies, yet the ubiquity of this compo- 
nent in other disks remained largely unrecognized until 
the UV observations. 

This neglect stemmed in part from the fact that 
Ha traces a limited subpopulation of the outer disk, 
namely those regions with O and B stars. Even 
relatively young regions may lack OB stars because 
they are older than 10 Myr or because they simply 
did not form such massive stars (Ferguson et al. "19981 
Herbert-Fort et al. 2009; Pflamm-Altcnburg & Kroupa] 
20081 ). As such, the ubiquity of outer disk star formation 
is somewhat concealed. These barriers are removed with 
UV observations and we now know that many nearby 
galax ies host young outer disk stellar populations (> 
30%; iThilker et alll2007l : IZaritskv fc Christlein I [200% . 
While galactic interactions can dramatically increase 
the level of outer disk star formation (see the well- 
known cases of M83, NGC 4625 and M94: IThilker et al. I 
[20051: IGil de Paz et al. l[200l iTruiillo et al. II2009D ." even 
isolated galaxies show low-level s of ongoing st a r for- 
mation in their outer parts (Ferguso n et al. 119981 : 
iChristlein fc Zarits kv 2008; Herbert-Fo rt et al. I [2009f ). 
However, the Ha knots provide the bright emission lines 
that make it possible to measure the kinematics we 
present here. 

While the gaseous outer dis k of NGC 628 has bee n 
well-studied, for example by Ka mphuis fc Briggsl (|1992| ). 
the stellar component near the outskirts of the disk re- 
mains poorly characterized. Observers have acquired 
deep broadband imaging of the outer disk of NGC 628 
(e.g. iNatali et al. I Il992i ). however the diffuse light be- 
comes difficult to reliably trace fainter than ^ 28 mag 
arcsec"^ in V, or beyond ~ 1.3i?25- To sidestep the dif- 
ficulty of obtaining detailed kinematic s from low surface 
brightness emission, [Herrmann et al. 1 (|2009l ) use plane- 
tary nebulae (PNe) across the face of NGC 628 to trace 
the underlying stellar distribution and to estimate the 
total mas s density of th e disk using the kinematic ap- 
proach of Ivan der Kruit [ (|1988f ) (specifically, by measur- 
ing the vertical velocity dispersion of a tracer of the 
disk mass distribution to estimate the underlying sur- 
face mass density). In principle, the PNe provide an 
excellent approach to the problem, with large numbers 
of PNe av ailable for analysi s . Unf ortunately for our pur- 
poses, the [Herrmann et al. I ([2009 ^ study of PNe in NGC 
628 is mostly constrained to the inner disk (only two of 
their PNe are beyond R25 ; their study of M83 and espe- 
cially M94 yield more PNe in those outer disks, however 
those are interacting systems for which the kinematics 
are complicated by recent events). The current work 
represents the first to analyze the disk kinematics out 
to 1.8i?2.'s in NGC 628 . Nevertheless, the results from 
[Herrmann et al. [ (|2009( ) at the largest radii R25) pro- 
vide an interesting comparison to the results found here, 
especially because the PNe trace a population of objects 
that is on average at least 100 x older than the Ha knots 
considered here. 

We use multiobject spectroscopy to search for Ha emis- 
sion associated with blue GALEX sources in the outer 
disk of NGC 628. We measure the dispersion of the Ha 
knot vertical velocity distribution, estimate the surface 
density of the outer disk, and determine if the baryonic 
material can account for the inferred mass. We then com- 



pare the results from the young Ha knots to those from 
the PNe and consider the implications on disk evolution 
from the two measurements. Section 2 presents our sam- 
ple selection, observations and data reductions. Section 
3 presents our analysis, results and discussion, including 
our measurement of the vertical velocity dispersion of the 
Ha knots and the surface mass density of the outer disk. 
Section 4 presents a summary and our conclusions. 

2. SAMPLE SELECTION, OBSERVATIONS, AND DATA 
REDUCTION 

We selected our target sample from sources 
identified in GALEX iin aging of NGC 628 by 
[Zaritskv fc ChristleiiTI (|2007[ ). We target 'blue' sources 
with FUV - NUV < 1 and NUV < 25, corresponding 
to the expected colors of star clusters with ages < 360 
Myr. We observed two fields that cover the outer disk 
of NGC 628 with the Inam ori Magellan Areal C amera 
and Spectrograph (IMACS: [Dressier et al~|[2006[ ) using 
the 'short' f/2 camera mode and the 600 line grating 
centered at 7695 A, providing a spectral resolution of 
~ 0.6 A pixel"^, or ~ 27 km s""'^ pixel^^ near Ha, on 
August 2, 3, and 4, 2005 (UTD). The IMACS field 
of view in the short camera mode is ^ 27' on a side, 
providing coverage out to 3i?25 {R25 = 5.2 arcmin; 
[Kamphuis fc Briggs I [1991 . The IMACS plate scale in 
our setup is 0.2 arcsec pixel"^. We cut 162 and 151 slits 
in the masks of the northeastern and southwestern fields, 
respectively, giving priority to the brightest GALEX 
sources with the bluest FUV — NUV colors (clusters 
with Ha emission are expected to have ages < 10 Myr, 
powered by massive stars that will make the clusters 
blue in the GALEX bands). The seeing remained below 
1" during the observations. We observed the northern 
field for a total of 18000 seconds, or 5 hours, and the 
southern field for half that time. To be certain that no 
significant mask/sky alignment drifts occurred during 
the overall integrations, we observed neither field for 
more than 2700 seconds in a single exposure. We 
observed HeNeAr comparison arc lamps through each 
mask for wavelength calibration. 

We reduce the data with the Carnegie Observato- 
ries System for MuhiObject Spectroscopy (COSMOS; A. 
Oemler et al.) package and extract individual 2D recti- 
fied, wavelength-calibrated and sky-subtracted spectra 
in the wavelength range 5600 - 9000 A. Using standard 
IRAF routines^, we median combine the individual sky- 
subtracted spectra corresponding to each knot, rejecting 
pixel values > 3(T from the median of the pixel stack to 
eliminate cosmic rays. 

3. ANALYSIS, RESULTS & DISCUSSION 

3.1. Ha detections 

We begin our spectral analysis by visually inspect- 
ing all of the sky-subtracted and combined spectra to 
search for detectable Ha flux from the targeted UV- 
bright knots. For any cases with visually-detected Ha 
emission, we average the rows showing signal and mea- 
sure the radial velocity centroid from this final ID spec- 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



trum. The procedure is designed such that we may ex- 
clude any rows that have minor overlaps with neighbor- 
ing spectra or contain any remaining artifacts. 

We estimate the line centroid measurement errors to be 
- 8 km s^i from fits to the bright OH 6498.72 A sky line 
from numerous (pre sky-subtracted) spectra. This mea- 
surement error is a lower limit to the Ha line measure- 
ment error because the sky line is typically brighter and 
cleaner than the Ha detections in our sample. We apply 
a constant —28.4 km s""'^ heliocentric velocity correction 
to the radial velocities of the Ha knots; the difference 
in the heliocentric correction between consecutive nights 
is ~ 0.1 km s~^, or ^ 1% of the typical line centroid 
measurement error. 

Of the 313 slits on the two IMACS masks, - 25% re- 
turned unusable spectra due to overlapping spectra or 
other artifacts in the data. We therefore only search for 
Ha emission in the remaining 235 spectra. Fifteen of 
these spectra (or ~ 6%) yield detectable Ha emission. 
Figure [1] presents a few of our final ID spectra illustrat- 
ing the range of strength of the Ha detections. Our 6% 
return is roughly what is expected given the difference in 
ages between the 'blue' GALEX-detected knots targeted 
here (< 360 Myr) and Ha knots 10 Myr) - especially 
when considering that priority was given to the bluest 
and brightest GALEX knots (so most targeted GALEX 
knots are likely < 200 Myr old). We plot the locations of 
these fifteen knots in Figure [5] and list the relevant infor- 
mation in Table [TJ The physical scale of the individual 
Ha knots (calculated from the extent of Ha emission 
along the slit) range from 25 to 92 pc, with a median 
value of 67 pc (adopting a distance of 8.6 Mpc to NGC 
628, dete rmined using the plan etary nebula luminosity 
function; [Herrmann et al.l 120081 ). 













670) 




















6570 6572 



6574 6576 6578 _ 6580 
Observed Wavelength (A) 



6582 



Fig. 1. — Three of our final ID spectra, illustrating the range 
of strength of the Ha detections. Note the different j/-axis ranges. 
QuaUtatively, our final sample is composed of two 'Strong', nine 
'Typical', and four 'Weak' Ha detections. 



3.2. Ha knot velocities in the disk frame and Uz 

Our present aim is to measure the vertical veloc- 
ity dispersion of the knots in the outer disk of NGC 




Fig. 2. — Figure 5 of IKamphuis fc Briggs I l|1992) showing the 



velocity field of the H I-emitting gas disk surrounding NGC 628, 
with our 15 Ha detections overplotted as encircled filled squares. 
The image is ~ 30' on a side. The Ha knots are found between 
~ 1 and 1.8 i?25 (13.5 - 23.2 kpc) 



TABLE 1 

Ha DETECTIONS 



RA (J2000) 


Dec (J2000) 


R/R25 


'^helio 


vth 


Ad 


01:36:04.47 


-1-15:45:27.3 


1.76 


640 


641 


-1 


01:36:05.40 


-1-15:44:23.3 


1.76 


668 


642 


26 


01:36:06.53 


-1-15:49:37.0 


1.72 


652 


647 


5 


01:36:14.07 


-1-15:48:19.0 


1.32 


667 


652 


15 


01:36:20.94 


-1-15:43:33.5 


1.17 


644: 


643 


1 


01:36:22.12 


-(-15:42:58.8 


1.20 


628 


643 


-15 


01:36:38.71 


-1-15:40:52.1 


1.19 


651 


642 


9 


01:36:51.02 


-1-15:52:03.6 


1.06 


680 


684 


-4 


01:36:59.23 


-1-15:52:13.6 


1.29 


657 


670 


-13 


01:37:00.35 


-1-15:51:19.1 


1.20 


659 


675 


-16 


01:37:00.94 


-1-15:49:49.7 


1.04 


658 


674 


-16 


01:37:05.30 


+15:44:45.3 


1.18 


639 


649 


-10 


01:37:09.28 


-1-15:41:14.0 


1.70 


704 


679 


25 


01:37:10.06 


-(-15:52:08.9 


1.64 


664 


677 


-13 


01:37:17.86 


-(-15:43:51.0 


1.79 


656 


711 


-55 



all velocities in km s~ . t>heiio refers to the heliocentric ve- 
locities of the Ha knots, while vth refers to THINGS heliocen - 
tric H I-emitting gas velocities measured bv lWalter et al.l II2008I ). 

628 (cr^, with z defined perpendicular to the disk 
plane). Given the low inc lination of NGC 628 (i ^ 6°; 
IKamphuis fc Briggs 11199^ , the observed line-of-sight ve- 
locities are almost ent irely due to the actua l vertical ve- 
locities. However, as IKamphuis fc Briggs I (|T992f) have 
shown, the outer gas disk of NGC 628 is chaotic and 
non-uniformly rotating. Whether this is also true for 
the outer stellar disk is unknown, and therefore whether 
we should use the previous orientation parameters ob- 
tained for the inner disk or some other ones is unclear. 
We therefore present results using various models for the 
disk orientation and kinematics, and accept the model 
with the lowest resulting velocity dispersion as the one 
most likely to be correct. Results of the following analy- 
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sis will be discussed further in §4. 

We use a maximum likelihood approach to estimate the 
Gaussian parameters best describing the unbinned rela- 
tive velocity (Aw) distributions. Specifically, we search 
a grid of mean and dispersion values (from —50 to 50 
km s^^ and 1 to 100 km s~^, respectively) for the com- 
bination that yields the highest total probability (likeli- 
hood) of the Gaussian function, using the unbinned Av 
values as the independent variable. Our velocity disper- 
sion estimates will provide only a single value for the 
radial range spanned by the knots (~ 1 — 1.8i?25)- We 
have insufficient detections to measure the velocity dis- 
persion as a function of radius. The following analysis 
will therefore consider the outer disk as a constant ve- 
locity dispersion, fixed-height component, and the values 
derived only represent the disk between ^ 1 — 1.8i?25- In 
support of this approach , note that the simulations of 
iKazantzidis et ahl (|2009[ ) predict a roughly constant ve- 
locity dispersion in the outer disk. 

We examine four models for the outer disk orientation 
and kinematics. First, we test a purely face-on disk by 
assuming a constant 656 km systemic velocity. This 
'model' yields a velocity dispersion of 18 ± 5 km s^^. 
Second, we adopt t he parameters of the i nner neutral 
hydrogen gas disk of iKamphuis &: Briggs I (|1992D (incli- 
nation i — 6.5°, major axis position angle PA = 25° 
and a mean velocity of 656 km s~]_]_and a rotation ve- 
locity of 170 km s~i (Fathi et al. Il2007f ) and find a ve- 
locity dispersion of 19 ± 5 km s~^. Third, we adopt 
the two-component model of lKamphuis fc Briggs I ()1992f ) 
(i = 13.5° and PA = 75° beyond 1.2i?25, using the same 
inner disk values as the previous case for the knots found 
at i? < 1.2i?25) and find a velocity dispersion of 31 ± 10 
km s~^. Finally, we adopt the neutral gas velocities from 
TheH I Nearby Galaxy Survey (THINGS; IWalter et al.1 
|2008| ) at the position of each knot and find a velocity dis- 
persion of 19 ± 6 km s^^. Figure [2] shows the spatial dis- 
tribution of the knots around the gas disk of NGC 628, 
and THINGS gas velocities at the positions of each knot 
are listed in Table [T] This 'model' tests for the possibil- 
ity that the knots have low velocity dispersion around a 
complex velocity field defined by the H I-emitting gas. 
Given the result that the velocity dispersion is similar 
to the simple disk models, we conclude that there is no 
evidence for highly complex kinematics. 

The velocity dispersion estimates generally cluster 
around 20 km s""'^ and agree to within 1 a (although 
the one estimate at 31 ± 10 km s^^ is just consistent 
with the others). The general agreement, independent of 
the model used, suggests that the dominant source of un- 
certainty is currently the statistical one associated with 
the limited sample size. Because the fourth 'model' is 
taken directly from the H I data and produces a velocity 
dispersion as tight as do the other models, we adopt this 
model for the remainder of the paper. 

Figu re [H shows th e distr ibution of velocities relative 
to the IWalter et al. I ()2008[ ) gas disk values (the fourth 
model), binned by 15 km s~^. We removed a slight, 
but significant, mean velocity offset (17 km s^^) when 
evaluating the velocity dispersion; Table [1] and Figure |3] 
include this correction. We attribute this shift to a veloc- 
ity zero point difference in the various datasets although 
we were unable to directly identify the cause (Figure [2] 



shows that the Ha knots are well-distributed across the 
face of the disk; we therefore do not expect that the off- 
set is due to an undersampling of the relative velocity 
distribution). When we exclude the knot at 01:37:17.86, 
-1-15:43:51.0 with Av = VheUo — vth = —55 km s~^, 
which lies on the high veloc ity cloud (HVC) described in 
IKamphuis fc Briggs I ()1992| ) , we arrive at our best (lowest 
and most likely) estimate of the observed Ha knot veloc- 
ity dispersion, 14 ± 4 km s~^. Excluding the same knot 
from the other models does not significantly change their 
values. The Gaussian representing this final unbinned Av 
distribution is overplotted in Figure [3] (for presentation, 
the height of the Gaussian is determined by a fit to the 
binned Av distribution; the dispersion is the only value 
used in our analysis, however). 
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Fig. 3. — Distributi on of Ha knot vertical velocities, relative 
to THINGS ( Walter e t al. 11200^ ') neutral gas measurements. The 
outlier is the source detected on the high velocity cloud; we ignore 
this source in the analysis because its relation to the gas disk is 
uncertain. 

We now consider the degree to which our observa- 
tional errors inflate the velocity dispersion measurement 
(cTofcs = 14kms~^). Subtracting our previous velocity er- 
ror estimate (8 km s~^) in quadrature from the measured 
velocity dispersion, we estimate that cr^ < 11 km s~^. We 
emphasize that is likely to be an upper limit to the 
vertical velocity dispersion because 1) the measurement 
uncertainty is larger than 8 km s^^, 2) the uncertainty 
in the reference H I-emitting gas velocities was ignored, 
and 3) we may have missed the correct model for the 
underlying disk kinematics. 

Our low velocity dispersion measurement suggests a 
relatively undistur bed stellar disk, which at fi rst seems 
to contrast with the lKamphuis fc Briggs ' ('1992') observa- 
tionof a chaotic outer gas disk. The Kamphuis fc Briggs] 
(1992) observation was strongly influenced by the HVCs 
and kinematic features at the largest radii, however-the 
bulk of the gas disk is in fact rather undisturbed (see 
Figured]). Nevertheless, a difference between a relatively 
undisturbed stellar disk and a disturbed gas disk might 
be expected, co nsidering the results of iThilker et al. I 
(i2007,) and iMoffett et al. I pOOa ). who find that Type 
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1 XUV disks (of which NGC 628 is a member) cover a 
wide range of galaxy masses, colors and morphologies, 
and that these disks may in fact be associated with mi- 
nor interactions or accretion. The disturbed structure 
observed in the gas may result from an interaction event 
(jKamphuis &: Briggs 1992) to which the stellar disk was 
less sensitive. Such an event may have then triggered the 
Type 1 XUV disk that is observed today. 

3.3. Outer disk scale height 

Assuming that the Ha knot density has the same ver- 
tical distribution as the underlying mass density distri- 
bution of the disk, for an adopted scaleheight we con- 
strain the surface mass density of the outer disk S using 
our measurement of cr^. We will return to discuss the 
validity of this assumption, although we do not have any 
alternative given our current data. The relevant equation 
(|van der Kruit II1988D is 

al = KG^K_, (1) 

where K = 1.705l7r for an 'intermediate' vertical mass 
density distribution (between the isothermal sech^(z) 
and exponential distributions) and G is the gravitational 
constant. The relation was constructed by solving the 
Boltzmann and Poisson equations for plane parallel lay- 
ers. For a galaxy with a flat rotation curve, the plane 
parallel case provides a realistic mo del from which to es - 
timate t he surface mass den sity (van der Kruitl[l988h . 
Because iChristlein fc Zaritskv (20081 ) found mostly flat 
rotation curves in a sample of 17 nearby outer disks, we 
expect the plane parallel framework to remain valid in 
the outer regions (we further assume a constant scale 
height over the radial range of the knots). Corrections 
to the re lation from a dark ha lo term are negligibly small 
at low z (|van der Kruit IIT9881 ) : however, in the outer disk 
the baryonic disk contribution to the mass budget pre- 
sumably drops more rapidly with radius than does the 
dark matter contribution, and so the dark matter may in 
fact not be negligible in t he outer disk. I n the following 
section we show that the Ivan der Kruit I (1988) relation 
is indeed appropriate at these radii (we find that no dark 
matter in the disk is required to explain existing obser- 
vations, even when assuming a rather large value for hz). 

Because NGC 628 is nearly face-on, we must estimate 
the scale hei ght hz of the t racer particles (the knots) 
indirectly. H errmann et al. I ()2009l ) argue, based on sta- 
bility arguments and the observed velocity dispersion of 
PNe, that the scaleheight near the edge of the region 
they probe 11 kpc) must be > 400 pc. The ver- 
tical disk flaring suggested by the PNe analysis resem- 
bles the known flaring of gas disks in their outer regions, 
with neutral hydrogen gas scale heights near ^ 700 pc at 
the R/ R^f, distances cons i dered here (IMerrifield 1 119921 : 
ICorbelh fc Saloeter I [1991 lOlling 1 119961). Because the 
Ha knots are young and because they lie at radii larger 
than those probed by the PNe, we adopt a scale height for 
the knots in line with that of the flaring gas, hz = 700 pc 
at 1.3i?25, the median radius of our Ha detections. Re- 
call however that we treat the outer disk as a fixed-height 
component. In reality, hz may range between 400 pc 
and ^ 1000 pc over the radii covered by our detections. 
Assuming a constant 700 pc over the entire outer disk 
therefore implies a large uncertainty in the final outer 



disk mass, which we assume to be ^ 20%. This level 
of uncertainty will not significantly affect our conclu- 
sions, however. Finally, we note that our assumption 
that the Ha knots trace the same mass distribution as 
the outer disk gas is different from what one would ex- 
pect in the inner disk. In the inner disk, molecular clouds 
(and their Ha knots) have a smaller vertical scaleheight 
than the mass, which is dominated by the stellar com- 
ponent. In the outer disk, however, where the dominant 
baryonic component is the gas (as we show below), our 
assumption is much more appropriate. Whether the Ha 
knots trace the neutral gas is more uncertain, although 
where examined in detai l the correspondence app ears 
good (jThilker et al. |[2005HHCTbert-Fort et al. II2009D . In 
view of this, a large scaleheight for the outer stellar disk 
should be expected. We caution however that if a sig- 
nificant amount of mass extends to heights larger than 
the 700 pc scaleheight assumed for the knots, we are not 
sensitive to it with the Ha knot measurement. 



3.4. Outer disk mass density 

We now estimate the mass density of the outer disk in 
order to determine the relative contributions of stars and 
gas to the total outer disk mass. Using CTz = 11 km s^^ 
and our adopted hz — 700 pc, we calculate that the outer 
disk surface mass density S is 7.5 M© pc~^. The total 
mass between 1 — 1.8i?25 is then ~ 8.4 x lO^Af©. 

We now use results from deep optical imaging of the 
oute r regions of NG C 628 {^v ^ 27 mag arcsec"'^ at 1.3 
i?25; iNatali et al. I [l992 ) to estimate that the mass-to- 
light ratio T of the outer disk is ~ 18. The outer disk 
mass is therefore not p rimarily in the fo r m of normal 
stars. For comparison, iHerrmann et al. 1 (|2009f ) found 
that the PNe velocity dispersions imply T = 1.4 for the 
inner disk. Reconciling the outer and inner disk T values 
would require a factor of 10 changes in cr^ or S, neither of 
which is that poorly constrained. Therefore, as one might 
expect, we detect a significantly higher T in the outer 
regions. If we adopt T = 3 a s an upper limit for an old 
outer disk stellar population ()Bell fc de Jong 1120011) . we 
find that the dark material in the outer disk must be at 
least five times more massive than the stellar component. 
In order to match the ~ 27 mag arcsec"^ observation, 
an outer disk stellar population with T < 3 would require 
< 1.3Mo pc-2 and < 1.4 x IO^Mq. This leaves 
> 7 X IO^A/q of dark material to be accounted for. 

We now consider the neutral gas component of the 
outer disk and whether it can account for the optically 
non-luminous mass. Using the out e r disk H I-emitting 
gas mass from iKamphuis fc Briggs! (|1992f ) (over the rel- 
evant radial range), ~ 6 x 10^ Mq, we estimate the to- 
tal outer disk baryonic mass (atomic gas -I- stars) to be 
between ^ 6 x 10^ — 7.4 x 10^ Mq. Comparing this to- 
tal baryonic mass to the total mass required by the Ha 
knot velocity dispersion (~ 8.4 x IO^Mq), we find that 
the baryonic matter (mainly the H I-emitting gas) can 
account for nearly all of the measured outer disk mass 
(note that we have ignored any molecular gas contribu- 
tion). Our assumption that the dark matter contribution 
to the disk is negligible at large radii appears justified, 
as does our claim that the gaseous component dominates 
the baryon budget. 
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3.5. Star formation history constraints 

A separate intriguing question is whether star forma- 
tion has occurred at the current rate for the entire hfe- 
time of this galaxy. The resuh of such a stellar com- 
ponent could, in principle, violate the surface brightness 
or surface mass density measurements. Independent of 
the surface brightness measurement we can simply cal- 
culate whether the integrated star formation results in 
a mass density in conflict with the measurement. Us- 
ing deep Ha imaging of three nearby late-type spiral 
galaxies (one of which is the current focu s ; NG C 628, 
NGC 1058 and NGC 6946'l. lFerguson et al. I (Il99l found 
their outer disk star formation rate densities to be be- 
tween 0.01 - O.OSM© pc"2 Gyr-^ For the follow- 
ing exercise, we adopt a SFR density of O.OSM© pc~^ 
Gy r~^ for the outer disk o f NGC 628. For comparison to 
thelFerguson et al. I (|1998[ ) SFR densities. PiVuiillo et al. I 
(12009D used GALEX FUV imaging to measure the SFR 
density in the outer disk of M94 and found ~ QAMq pc~^ 
Gyr~^. The star formation in the outer regions of M94 
is unusually large however, producing an anti-truncated 
optical disk light profile (NGC 628 shows a standard ex- 
ponential optical light profile to i?25, with only a very 
slight possible downbending in the profile to ~ 1.3i?2.'^, 
the extent of the deep optical observations: iNatali et alH 
[19921 For NGC 628, assuming uniformly-distributed 
constant star formation over a Hubble time leads to an 
expected stellar outer disk mass of 5 x lO^M© (0.45Mo 
pc~^), well within our measured mass of 8.4 x lO^M© 
and within the stellar mass inferred from the measured 
surface brightness (M* < 1.4 x lO^M© for a population 
with T < 3). We conclude that the outer disk may have 
been forming stars for the entire lifetime of the galaxy at 
approximately the current star formation rate. 

3.6. Similarities between the young and old components 

We now compare the outer disk velocity dispersions 
of the young (< 10 Myr) and old (100 Myr - 10 
Gyr) populations, trac ed by the Ha knots and PNe 
(jHerrmann et al.ll2009[ ). respectively, and find little dif- 
ference (~ 11 km s""'^ vs. ~ 12 km s~^; a small difference 
may be allowed if we treat the Ha velocity dispersion 
as a true upper limit, however). Because the difference 
is unlikely to be much larger than the la error of the 
measurement, we cannot say that there is a significant 
difference in the kinematic or structural properties be- 
tween the young and old outer disk populations. 

A powerful tracer of the evolution of stellar populations 
of different ages is the increase of velocity dispersion with 
age, or, alternatively, the varying asymmetric drift of var- 
ious populations. The observed low az of the Ha knots 
immediately suggests that the degree of asymmetric drift 
will be small. If one considers that some fraction of the 
observed value must reflect the intrinsic dispersion of the 
gas, then the asymmetric drift will be even smaller. The 
other population that has been measured, PNe, also has a 
low velocity dispersion and the inferred asymmetric drift 
is als o relatively small (< 15 km s~^; [Herrmann et al. 1 
I2009D . With the current measurements we cannot con- 
clude that there has been any increase in the velocity 
dispersions between the very young Ha knots and the 
older PNe (we comment further on this issue in the fol- 
lowing section). It is evident that such an investigation 



will require uncertainties in the velocity dispersions bet- 
ter than ^ 1 km s^^. 

4. SUMMARY AND CONCLUSIONS 

We have measured the kinematics of Ha knots in the 
outer disk of NGC 628. We find the stellar disk (traced 
by the Ha knots) to have a low velocity dispersion, 
suggesting an undisturbed extended stellar disk. NGC 
628 differs from the better-known outer disks of M83 
and M94, which were likely accentuated by recent dis- 
turbances from neighbors, but resembles (both in op- 
tical and kinematic profiles) the outer disks detected 
in other relatively und i sturbe d nearby edge-on disks by 
iChristlein fc Zaritskv I ()2008D . supporting the idea that 
outer disk star formation can be a low-level and ongoing 
phenomenon in isolated galaxies. 

We find (Tz of the Ha knots to be < 11 km s~^ be- 
tween 1 — 1.8i?25 (13.5 — 23.2 kpc). We adopt a scale- 
height similar to the known flaring gas profiles of outer 
disks {hz = 700 pc) and estimate a mass density E = 7.5 
M0 pc~^ that can be entirely explained by the observed 
gas and stars in the outer disk. If the Ha-hosting disk 
is actually much thinner than the flaring gas disk, more 
dark matter in the outer regions would be allowed. As- 
suming that outer disk star formation has been going for 
a Hubble time does not violate either the surface bright- 
ness nor surface mass const raints. The high incidence 
of outer disks (c.f. Christ lein fc Zaritskv 1 12008') suggests 
that the star formation is not a rare phenomenon - here 
we show that current limits cannot exclude long-lived 
outer disk star formation. 

Finally, the velocity dispersio n of PN e towards the 
outer disk of NGC 628 (Herrman n" et al. | |20091 is nearly 
the same as that of the Ha knots (the discrepancy grows 
slightly when considering the dispersion of the Ha knots 
as a strict upper limit, although any plausible difference 
remains small). This can result either if the PN popu- 
lation is rather young (so that scattering has not had a 
chance to enlarge the dispersion significantly), if there 
is very little scattering, or if scattering occurs primarily 
in a single (or few) discrete events that occurred prior 
to the creation of the bulk of the PNe. We do not ex- 
pect the outer disk to be exclusive young (see above) and 
we do expect some level of scattering in outer disks (if 
not from the classical spiral arms and molecular clouds 
of inner disks, then perhaps from satellites and dark 
halo substructure) . The solution may be found in the 
iKazantzidis et ahl (|2009f) simulations, which show that 
outer disk heating is dominated by the most massive in- 
fall event of halo substructure onto the disk, so that stel- 
lar populations of different ages do not necessarily have 
different velocity dispersions. A larger sample of outer 
disk kinematic measurements could be used to constrain 
the rate and impact of such infall events. 
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